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In this article, the prospects for studying the tensor structure of the HZZ vertex with the LHC 
experiments are presented. The structure of tensor couplings in Higgs di-boson decays is investigated 
by measuring the asymmetries and by striding the shapes of the final state angular distributions. 

The expected background contributions, detector resolution, and trigger and selection efficiencies 
are taken into account. The potential of the LHC experiments to discover sizeable non-Standard 
Model contributions to the HZZ vertex with 300 flV 1 and 3000 fb” 1 is demonstrated. 


I. INTRODUCTION 

In the Summer of 2012, the CMS and ATLAS Collab¬ 
orations at the LHC reported the discovery of a new neu¬ 
tral resonance in searches for the Standard Model Higgs 
boson. This discovery was later confirmed by analyses 
of the full LHC Run-I dataset by both collaborations 
[1, 2]. It was demonstrated that the new particle with 
a mass around 125.5 GeV was dominantly produced via 
the gluon-fusion process and decays into pairs of gauge 
bosons: 77 , ZZ and WW. The observed production and 
decay modes identified the discovered particle as a neu¬ 
tral boson. The subsequent measurement of its couplings 
to fermions and bosons demonstrated the compatibility 
of the discovered resonance with the expectations for the 
Standard Model Higgs boson within available statistics 

[3-5]- 

In the Standard Model, electroweak symmetry break¬ 
ing via the Higgs mechanism requires the presence of a 
single neutral Higgs boson with spin 0 and even CP- 
parity. Theories beyond the Standard Model often re¬ 
quire an extended Higgs sector featuring several neutral 
Higgs bosons of both even and odd CP-parity. In such 
a case, mixing between Higgs boson CP-eigenstates is 
possible. The Higgs boson mass eigenstates observed in 
experiment may thus have mixed CP-parity. Such an ex¬ 
tension of the Higgs sector is important because effects of 
CP violation in the SM are too small and, in particular, 
cannot explain the baryon asymmetry of the Universe. 
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Dedicated studies of spin and parity of the Higgs can¬ 
didate discovered by ATLAS and CMS showed that its 
dominant spin and parity are compatible with J CP = 
0 ++ [4-6]. The dataset of about 25 fb -1 currently col¬ 
lected by each of the major LHC experiments allows 
to set an upper limit on the possible CP-odd contribu¬ 
tion. The sensitivity is expected to improve with larger 
datasets to be collected at the LHC. 

There have been many works on direct measurement 
of CP violation in the Higgs sector [7 -32]. In this paper 
the sensitivity of LHC experiments to observe CP-mixing 
effects with 300 fb -1 and 3000 RU 1 is evaluated using the 
method of angular asymmetries. 

This paper is organised as follows. In Section II ob¬ 
servables sensitive to CP violation in the HZZ vertex are 
discussed. The spin-0 model, a Monte Carlo production 
of signal and background, and a lagrangian parametrisa- 
tion for CP-mixing measurements are discussed in Sec¬ 
tion III. In Section IV the expected sensitivity of the LHC 
experiments to the CP-violation effects based on angu¬ 
lar asymmetries is presented. Constraits are set on the 
contribution of anomalous couplings to the HZZ vertex. 
Section V introduces the measurement technique based 
on observables fit. Exclusion regions for the mixing angle 
are presented. Section VI gives the overall summary of 
obtained results. 

II. OBSERVABLES 

In this paper we study the sensitivity of final state ob¬ 
servables to the CP violating HZZ vertex in the process: 

gg -A H ZZ -> 41. (1) 

Following the notation introduced in [21], the general 
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scattering amplitude describing interactions of a spin- 
zero boson with the gauge bosons is given by: 

A(X -7 VV) = l{gim 2 v ete* 2 +g 2 f;^r^ 

+ gif$ ) r WlM '). ( 2 ) 

Here the f^» lv = — e/gf is the field strength tensor 

of a gauge boson with momentum g, and polarisation 
vector €i] = l/ 2 e /il/a ^/ a/ 3 is the conjugate field 

strength tensor. The symbols v and my denote the SM 
vacuum expectation value of the Higgs field and the mass 
of the gauge boson respectively. 

In the Standard Model, the only non-vanishing cou¬ 
pling of the Higgs to ZZ or WW boson pairs at tree- 
level is gi = 2 i, while g 2 is generated through radiative 
corrections. For final states with at least one massless 
gauge boson, such as 77 , gg or Z'y, the SM interactions 
with the Higgs boson are loop-induced. These interac¬ 
tions are described by the coupling g 2 . The coupling 54 
is associated with the interaction of CP-odd Higgs boson 
with a pair of gauge bosons. The simultaneous presence 
of CP-even terms g\ and/or g 2 and the CP-odd term p 4 
leads to CP violation. 

In general, gt couplings can be complex and momen¬ 
tum dependent. However imaginary parts of these cou¬ 
plings are generated by absorptive parts of the corre¬ 
sponding diagrams and expected to be small: approxi¬ 
mately less than 1%. We further assume that the energy 
scale of new physics is around 1 TeV or higher, so that 
the momentum dependence of the couplings can be ne¬ 
glected. Thus, in the following we will consider gi cou¬ 
plings as real and momentum-independent. 

These assumptions are entirely consistent with the 
framework of an effective field theory (EFT) of the SM. 
If the energy scale of the new physics is much higher than 
the electroweak scale new effects can be described by an 
EFT with the SM Lagrangian supplemented by higher 
dimension operators of d = 6 . Such an approach was 
worked out in detail in [33, 34] 

One of possibilities to study CP violation in the process 
of Eq. 1 is to analyse the shapes of angular and mass 
distributions of the final state [35, 36]. The common 
choice of angular observables for this type of analysis is 
show in Fig. 1. 

A complimentary approach is based on studies of 
angular-function asymmetries arising in the case of CP 
violation. There are six observable functions proposed in 
[37]. The first angular observable function is defined as 
follows: 


Q = ( V2Z -Piz) ■ {P3H + P4h) 

\P 2 Z — piz 11 P3H + P4H | 

Here p), i = 1,... 4 are the 3-nromenta of the final state 
leptons in the order l\l\l 2 l 2 . The subscripts Z and H 
denote that the corresponding 3-vector is taken in the Z 
or in the Higgs boson rest frames. Using these definitions, 



FIG. 1. Definitions of the CP-sensitive final state angular 
observables in gg —> H —> ZZ —>■ 4Z decay. 


the second observable function reads: 

0 _ (P 2 Z ~ Piz) • (P4H X p 3H ) 

|P 2 Z — piz | \P4H X P3H | 

The third observable function 0 3 is constructed using 
Or- 

O3 = Oi o 3a o 3b , 


where 

q _ (P4Z ~ P3Z) • (PlH X p 2H ) 

\P4Z — P3z\\PlH X P 2 h\ 

and 

_ (piz - P4z) ■ (pig +P 2 h) 

\P3Z - P4z\\piH + P 2 h\ 

The remaining three observable functions are given by: 

Q _ [{P3H X P4 H ) ■ Pih][(P3H X p iH ) ■ (pig X P2h)\ 

\P3H +P4H\ 2 \P1H +P2H\\P3Z ~ P4z\ 2 \PlZ -P2z| 2 /16 ’ 


0 5 


[( P4H X p 3H ) ■ pjg] [(pjz ~ P 2 Z) ■ P 3 Z] 
\~P3H + P4H\\P3Z - P4z\ 2 \P1Z — P 2 z | 2 /8 


and 

Q _ [(piz - P 2 z) ■ (P 3 H + P4h)][(P3H X p 4g ) • p 1H ] 
\piz - P 2 Z | : 2 | P3H + P4H | 2 |P 3 Z - P4Z | /4 


These observables are related to the final state angular 
variables defined in [35] and illustrated in Fig. 1. For 
instance, a trivial calculation yeilds: 0\ = cos and 
0 2 = — sin (j> sin Q\. 

Note that the total cross section is CP even (no inter¬ 
ference between CP-even and CP-odd terms) and cannot 
be used to detect the presence of CP violating terms in 
the HZZ vertex. 
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TABLE I. Definitions of effective tensor couplings gxvv' in¬ 
troduced in Eq. (3) in units of the Higgs vacuum expectation 

v. The symbol C is defined as: C = \J 


III. SPIN-0 MODEL AND MONTE CARLO 
PRODUCTION 

The dominant Higgs boson production mechanism at 
the LHC is gluon-fusion. To simulate the production 
of a Higgs-like boson and its consequent decay into ZZ 
and 4 1, the MadGraph5 Monte Carlo generator [38] was 
used. This generator implements the Higgs Characterisa¬ 
tion model [39]. The corresponding effective Lagrangian 
describing the interaction of the spin-0 Higgs-like boson 
with vector bosons is given by: 

^ = jc Q K S M [ ^9hzz Z^Z^ + gnww W+ W~ p ~\ 

^ \pO(K'H'y'i9H'y'y A^A^ T A^A^ J 

2 [pa^HZ'ydHZ'y Z^A^ T S^KAz^dAZ^Z^A^ 1 J 

- J [c a K Hgg g Hgg G%G a ^ + s a K Agg9Agg G“ u G a ^] 

^ 4 A [ CaKHZZ + s a n A zz Z^Z^] 

~ 2 A \. CaKHWW + s a n A ww W^W-^] 

- -^c a [khOj Z^dfj.A^ 1 ' + khoz Z^d^Z^ 

+ (n H dwW+d lx W-^ + h.c.)\ J>A, (3) 


does not assume that the Higgs boson belongs to a dou¬ 
blet of the weak SU( 2) group. Interaction terms corre¬ 
sponding to a Lagrangian of this type do not necessarily 
form a complete basis. However, this form is convenient 
for analysis of experimental data, as it relates in a simple 
way effective couplings and quantities observed in exper¬ 
iments. Note that there is a different and very popular 
EFT approach [34] to studies of the Higgs boson sector 
based on a complete set of operators of dimension six 
(the so-called Warsaw basis). 

The relations between parameters of the Lagrangian of 
Eq. (3) and tensor couplings of the effective amplitude of 
Eq. (2) can be derived from Feynman rules. The corre¬ 
sponding conversion coefficients are shown in Table II. 


Coupling 

ZZ 

WW 

77 

z 7 

gg 

gi/ZiCa 

ksM 

ksM 

- 

- 

- 

gil^iCa 

Ahzz 

Khww 

-^Of'7'7 

Ah z-y 

A’Hgg 

gi/2is a 

Kazz 

Kaww 

K-A^'y 

Kaz-i 

K A gg 

Si /2 iCa 

Khsz 

Re(KHdw) 

- 

- 

~ 

g’l /2ic a 

- 

ilm(KHdw) 

- 

- 

~ 


TABLE II. Conversion coefficients between parameters of the 
Lagrangian of Eq. (3) and tensor couplings of the effective 
amplitude of Eq. (2). 


In this table the following definitions are used: 

lv ~ , 

iAXVV’ ~ W^XVV KXVV i 


1 V ( Ai Y 


Kudv = oT - knov, 

2 A V my J 


where A is the new physics energy scale and the field 
strength tensors are defined as follows: 

Vp V = YK - d v v„ (F = A, Z, fF±), 

G a s~ia o . rabcs-ib s~ic 


c a = cos a, and s a = sin a. 

Here X denotes either H or A and the index VV denotes 
the final state gauge boson pair. The effective couplings 
gxvv’ are defined as follows: 

• In the case of ZZ or WW interactions, gxvv = 1; 


The dual tensor V jlv is defined as: 


V„ v = -^p°V p °. 


The mixing angle a allows the production and decay of 
CP-mixed states and implies CP violation when a / 0 
or a / 7 t/ 2. The definitions of effective tensor couplings 
gxvv are shown in Table I. 

The Lagrangian in Eq. (3) is an effective Lagrangian 
with U(1)em symmetry. It parametrizes all possible 
Lorentz structures, is not SU(2) x 17(1) invariant and 


• For 77 , Z^j and gg interactions, couplings gxvv 
are equivalent to the couplings gxvv defined in 
Table I. 

The couplings Khov , where V = IF, Z, 7 , correspond to 
the so-called contact terms of the Higgs Characterisation 
Lagrangian of Eq. (3). These contact terms can be re¬ 
produced in the amplitude of Eq. (2) by re-parametrising 
the gi coupling in the following form [40]: 

, (A ,2) n SM . '2 1 ^ 1 + 1^1 , 

9l\Ql^Q2)— 01 +9 1-T 2 -r 9i - ~j~2 -■ 
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This equation represents the leading terms of the form 
factor expansion. In the case of complex know, the mo¬ 
menta of the W bosons should be assigned as follows: 
qi for W~ and q 2 for W + . In the case of HZj interac¬ 
tion with a real photon, the term proportional to knoy 
vanishes. 

In the following we will consider a model based on the 
Lagrangian of Eq. (3) in which the mixing is provided 
by the simultaneous presence of the Standard model 
CP-even term and a non-Standard model CP-odd term 
in the HZZ decay vertex. The signal Monte Carlo 
samples used in this analysis are produced using the 
Higgs Characterisation model parameters presented in 
Table III. 


ksM 

kHZZ 

kAZZ 

^Hgg 

kAgg 

A, GeV 

1 

0 

28.6 

1 

1 

l(p 


TABLE III. Parameters of Higgs Characterisation model used 
for Monte-Carlo simulation of signal samples. 


The coefficient kAzz was chosen such that it pro¬ 
vided equal cross sections for decays of CP-odd and CP- 
even Higgs states: a(c a = 0) = a(c a = 1). The ten¬ 
sor couplings for the decay vertex corresponding to the 
amplitude of Eq. (2) can be restored using the follow¬ 
ing relations: g 2 = 2 ic a and g 4 = 2 is a KAzz, where 
Kazz = 1-76. It is noted that the factor 2 i is not impor¬ 
tant in the study of asymmetries because it defines the 
overall cross-section normalisation. 

The signal samples were produced using the Mad- 
Graph5 Monte Carlo generator [38]. These samples were 
created in the range of mixing angles — 1 < cos a < 1 
in steps of 0.05. The dominant background processes 
qq —> ZZ, Z 7 were also simulated with MadGraph5. 

After simulation of signal and background events at 
y/s = 14 TeV, the parton showering was performed us¬ 
ing the PYTHIA 6 Monte Carlo generator [41]. Generic 
detector effects were included by using the PCS pack¬ 
age [38]. The main detector parameters used for this 
simulation are presented in Table IV. For comparison, 


Parameter 

Value 

Electromagnetic calorimeter resolution 

■y/E 0.1 

Hadronic calolrimeter resolution- y/E 

0.8 

MET resolution 

0.2 

Outer radius of tracker (m) 

1.0 

Magnetic field (T) 

2.0 

Track finding efficiency 

0.98 

Tracking 77 coverage 

2.5 

e /7 7] coverage 

2.8 

Muon rj coverage 

2.8 


TABLE IV. Tuning parameters used to simulate detector ef¬ 
fects with PGS package. 

the expected acceptance, efficiencies and resolutions of 
the ATLAS and CMS detectors of the LHC can be found 
in [42, 43], 


Finally a kinematic selection was applied. It was required 
that candidates decayed to two same flavour oppositely 
charged lepton pairs. If several of such candidates could 
be reconstructed in an event, the leptons pairs with in¬ 
variant masses closest to the on-shell Z mass where cho¬ 
sen. Each individual lepton had a psudorapidity \r]\ < 2.5 
and transverse momentum pt > 7 GeV. The most ener¬ 
getic lepton should satisfy pt > 20 GeV whereas the sec¬ 
ond (third) similarly had pt > 15 GeV (pt > 10 GeV). 
The invariant mass of the on-shell Z boson was in the 
mass window (50,106) GeV while the off-shell Z boson 
inz* > 20 GeV. Only Higgs candidates in the signal re¬ 
gion 115 GeV< run < 130 GeV where considered. The 
selection is a simplified version of the one presented in 
[ 2 ]- 

IV. ASYMMETRIES 

For each observable O,; sensitive to CP violation, the 
corresponding asymmetry can be defined as: 

N{O t > 0) - N(Oi < 0) 

1 N(O t > 0) + N(Oi < 0) ’ 1 1 

where N is the number of events with the observable less 
or greater than zero. Integrating the corresponding de¬ 
cay probabilities, it can be shown that these asymmetries 
directly probe the tensor couplings defined in the ampli¬ 
tude of Eq. (2) [37]. The value of A 4 is proportional to 
Im{g 4 ), while A 2 , A 3 , A 4 , A 5 and A 6 probe the values of 
Re(g 4 ) and Im(g 2 ) respectively. 

Analysis of asymmetries sensitive to CP-violation for 
the process of Eq. (1) was performed in [37]. In this 
section we extend this analysis by including effects of 
parton showering, hadronization, generic detector effects 
and contributions from the irreducible qq —» ZZjZ 7 —> 
Al background. Lepton interference in the final state and 
the contribution of two off-shell Z-bosons are also taken 
into account. 

The distributions of observables 0 2 ,0 3 , O 4 and O 5 for 
two values of the mixing angle cos a = 1 and cos a = 0.5 
are shown in Fig. 2. Signal H —> ZZ —> Al events are 
generated using the production and decay model defined 
in Table III. The contributions from the signal and qq — > 
ZZ — > Al background are normalised to their respective 
expectations at 300/6“ V It is noted that the presence of 
CP-mixing leads to distortions of distributions of selected 
observables. The distributions of 0 2 through O 5 become 
asymmetric in the presence of a real component of g 4 . 
This asymmetry is especially pronounced for 0 4 and O 5 . 
As suggested in [37], the background is CP conserving 
and the corresponding distributions of observables are 
symmetric. The shapes of asymmetries A ; for the model 
presented in Table III are shown in Fig. 3. The pure 
CP-even and CP-odd cases are given by cos a = 1 and 
cos a = 0 , respectively. 

Note, that according to the structure of Lagrangian 
(Eq. (3)) the CP-violating contribution is defined by the 
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FIG. 2. Distributions of observables 02 , 03,04 and O 5 for 
two values of the mixing angle a. 



FIG. 3. Asymmetries generated for observables Oi. 


where N = Ng + Nb is the total number of signal and 
background events and AN is the difference in the num¬ 
ber of events with Oi < 0 and Oi > 0. It is also noted 
that AN ss A Ns, because the ZZ background does not 
contribute to asymmetries at leading order. Following 
the results of the simulation presented in [44], the num¬ 
ber of signal and background events at yfs = 14 TeV 
can be estimated as: Ng = 1.32Z, and Nb = 0.71 L 
respectively. Here L represents the integrated luminos¬ 
ity in fb _1 . A dataset with the integrated luminosity of 
300 fb ” 1 is expected to be collected during the Run III 
of the LHC. 

Using the above expressions, one can estimate an ex¬ 
pected asymmetry of about 9.5% to be measured with 
this data sample. The corresponding significance will 
be around two standard deviations. The region 0.340 < 
cos a < 0.789 will then be excluded at 95% CL. 

This exclusion range can be expressed in terms of f Qi 
fraction of events [4] arising from the anomalous coupling 
54: 


, = cdM 2 
l9i <u| 5 l | 2 +a 4 M 2 ’ 


(5) 


where gi are couplings of the decay vertex, and (J t is the 
cross section of the processes H —► ZZ —» 41 correspond¬ 
ing to gi = l,gi^tj = 0. Eq. (5) can be rewritten in terms 
of the mixing angle a as: 



parameter p = Kazz tana. This parameter thus deter¬ 
mines the corresponding asymmetries of angular observ¬ 
ables. Knowing the distribution of asymmetries for given 
Kazz it is possible to obtain the corresponding distribu¬ 
tions for any Kazz by using the condition p = const. 

It is noted, that for the physics model used in this 
study, the observables Oi and Oq do not generate asym¬ 
metries visible with the current Monte Carlo sample. 
The consistency of these asymmetries with zero confirms 
that additional effects that are taken into account in our 
work such as lepton interference, off-shell ZZ produc¬ 
tion, ZZ background, experimental cuts and detector 
acceptance do not produce an artificial asymmetry not 
related with the presence of CP-odd terms. The asym¬ 
metric behaviour is clearly visible for O 2 through O 5 . 
The asymmetries for O 4 and O 5 calculated using Eq. (4) 
may exceed 10 % . 

In Fig. 3 asymmetry plots are given for cos a in the 
range from 0 to 1. For negative cos a the asymmetries 
change sign but keep the same shape. This property al¬ 
lows using the asymmetry approach to measure the rela¬ 
tive phase in the amplitude of Eq. (2). 

The significance of the expected asymmetry can be es¬ 
timated as: 

S = AN/VN = AiNg/VN, 


where the ratio of cross sections 174/171 = 0.139 is ob¬ 
tained from the Monte Carlo generator. 

The range of the fraction of events of Eq. (5) close to 1 
has been already excluded by CMS [4]. Taking this into 
account, the exclusion limit obtained in the presented 
analysis becomes f 9i < 0.206 at 300 fb ” 1 for the model 
described by the Lagrangian of Eq. (3) and parameters 
given in Table III. 

For the high luminosity LHC, assuming the same signal 
and background yields per fb as above, the following ex¬ 
clusion range can be established: 0.089 < cos a < 0.968 
at 95% CL. This corresponds to an upper limit f 9i = 
0.028 at 3000 fb” 1 . 

In the same way as above, we performed estimates for 
four more values of the model parameter Kazz- Monte 
Carlo samples were generated for each point of two di¬ 
mensional model space (cos a, Kazz)- The number of 
signal events was calculated as Ng = Ng M a/a SM as¬ 
suming constant K-factors. The results are presented in 
Table V. These limits on / 94 are close to the ones ex¬ 
pected in ATLAS [44] and CMS [4] experiments. 

The region of -Kazz/ 1-76 above 1.4 is not considered. 
In this region the cross sections exceed the SM cross sec¬ 
tion by more than a factor of two. 

In Figs. 4 and 5 the regions of model parameter space 
(c Q , Kazz ) excluded by the current analysis are shown. 
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L, fb" 1 

300 

3000 

Kazz / 1-76 

^ Col 

f 94 

Z\Cq; 

fgi 

0.6 

- 

- 

0.122-0.921 

0.026 

0.8 

0.431-0.650 

0.274 

0.100-0.953 

0.027 

1.0 

0.340-0.789 

0.207 

0.089-0.968 

0.028 

1.2 

0.307-0.852 

0.191 

0.087-0.975 

0.031 

1.4 

0.297-0.886 

0.188 

0.086-0.981 

0.032 


TABLE V. Upper limit on f gi and cos a range excluded at 
the 95% CL. 

The shadowed areas are excluded at the 95% CL. Lines 
in Figs. 4 and 5 represent a polinomial fit to the results 
of the method of asymmetries. 

Note that CP-odd observables were studied also in [45]. 
According to this article the detection of CP-violating 
effects is out of reach of the LHC. However, as was men¬ 
tioned in [45], these effects might in principle attain large 
values because of numerical enhancements. 



FIG. 4. The 95% CL exclusion limits for model parameters 
c Q , Kazz at 300/6 -1 . Regions of rejected model parameters 
are shown. 



FIG. 5. The 95% CL exclusion limits for model parameters 
c Q , Kazz at 3000/fe -1 . Regions of rejected model parameters 
are shown. 


V. MIXING ANGLE OBSERVABLE FIT 

The asymmetries Ai discussed in the Section IV are 
integrated quantities of angular observables Oi and thus 
provide limited information about the anomalous con¬ 
tributions to the HZZ vertex. The optimal sensitivity 
to these contributions can be obtained by studying the 
shapes of distributions of observables Oi and their corre¬ 
lations. 

The sensitivity of individual observables to the pres¬ 
ence of anomalous contributions to the HZZ vertex is 
studied by fitting the shape of these observables as a 
function of the mixing angle. The likelihood function of 
the fit is defined as: 

Nchan Vi,i 

£(cosa,/z,0) = n n P(N i j\n j -S i j(cosa,8)+B l ' ] (e)). 

3 * 

Here besides the parameter of interest cos a, two nui¬ 
sance parameters have been introduced: the best fitting 
signal strength y, and a systematic normalization uncer¬ 
tainty 9. The likelihood function is a product over the 
different final states and bins of the specific observable 
that is being fitted. In each bin, the observed number of 
events from pseudo-data N, is compared to the expected 
number of events of the model S + B assuming a Pois- 
sonian distribution of entries P. By varying the mixing 
parameter cos a of the likelihood for a given dataset we 
can construct the standard log-likelihood test statistic: 

„. . . . „, £(cos a) 

—21nA(cosa) = —2 In —-—, 

£(cosa) 

where d denotes the mixing angle that maximises the 
likelihood function over the scan. The other likelihood 
parameters are profiled at the corresponding cos a value. 
The 95% exclusion is reached when —2 In A(cos a) > 3.84. 
The definitions of the 64% CL and 95% CL exclusion 
regions is demonstrated in Fig. 6. 

Results of the scan of the mixing angle a produced 
with the mixing angle observable fit corresponding to the 
integrated luminosity of 300 fb -1 are presented in Fig. 7. 
The results are reported for the model with Kazz = 1-76 
and remaining parameters as defined in Table III. The 
values of the mixing angle cos a used to generate the 
input pseudo-data are marked on the x-axis. Every bin 
of the injected cos a on represents the null-hypothesis 
likelihood curve similar to Fig. 6. The y- axis shows the 
cosd values reconstructed in the fit. The blue and grey 
dashed areas represent the 64% CL and 95% CL limits 
respectively. The white area in each bin of injected cos a 
is excluded at 95% CL. As expected, the sensitivity to 
the mixing angle varies for different observables, resulting 
in significantly different exclusion regions. The weakest 
exclusion is reached with the O 2 , while the strongest is 
reached with the O4. 

The results corresponding to the integrated luminos¬ 
ity of 3000 fb” 1 are presented in Fig. 8. Compared to 
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FIG. 6 . Example of the likelihood curve for the mixing an¬ 
gle observable fit of O 4 . The definitions of the 64% CL and 
95% CL exclusion regions are demonstrated. 


300 fb - , the 95% CL exclusion regions around the fitted 
cos a values are significantly reduced. Assuming the pure 
Standard Model signal, the following exclusion limits can 
be set using the O 4 observable alone: 0 < cosd < 0.708 
at the 95% CL for 300 fb -1 and 0 < cosd < 0.908 at 
the 95% CL for 3000 fb -1 . The exclusion limits obtained 
from other observables assuming the Standard Model sig¬ 
nal are reported in Table VI. 

The exclusion limits obtained for hypothetical BSM 
signals can be read from Fig. 7 and 8 . It is noted that by 
fitting the shape of the O 4 observable alone the exclusion 
limits similar to those reported in the Section IV can 
be obtained. Further improvements can be obtained by 
combining several observables in the same fit. 


L, fb -1 

300 

3000 

Observable 

A c a 

fg 4 

Z\Cq: 

fg 4 

Oi 

0 - 0.695 

0.315 

0 - 0.903 

0.089 

O 2 

- 

- 

0 - 0.604 

0.428 

Oa 

0 - 0.719 

0.287 

0 - 0.911 

0.081 

Or 

0 - 0.708 

0.300 

0 - 0.908 

0.084 

Os 

0 - 0.631 

0.394 

0 - 0.883 

0.108 

Oe 

0 - 0.533 

0.520 

0 - 0.852 

0.104 


TABLE VI. Upper limit on / 94 and cos a range excluded 
at the 95% CL with the mixing angle observable fit. The 
Standard Model signal is assumed. The BSM templates are 
generated according to the model defined in Table III with 
Kazz = 1.76. 


VI. CONCLUSION 

In this article studies of tensor structure of the HZZ 
vertex are presented. The investigation is performed us¬ 
ing the pp —> H — > ZZ —> 4 1 process assuming the gluon 
fusion production of the spin-0 resonance. The back¬ 
ground contributions, detector resolution, trigger and se¬ 
lection efficiencies expected for the LHC are taken into 
account. Two different approaches to detect CP-violation 
effects in the HZZ vertex were used. The first approach 
is based on a simple counting experiment for angular 
asymmetries of CP-sensitive observables. It was shown 
that the presence of CP violating terms may result in 
angular asymmetries exceeding 10%. The 95% CL exclu¬ 
sion ranges for the mixing angle at different parameters of 
spin-0 Higgs boson model including the Standard Model 
CP-even term and anomalous CP-odd term 34 are cal¬ 
culated. These results are also presented in terms of the 
effective cross section fraction f 94 . The obtained limits 
are comparable with the ATLAS and CMS projections 
for Run III at the LHC and the High-Luminosity LHC 
presented in [4, 44]. 

The sensitivity of individual observables to the pres¬ 
ence of anomalous contributions to the HZZ vertex was 
studied by fitting the shape of these observables as a 
function of the mixing angle. It is demonstrated that 
using a single most sensitive observable, this approach 
gives fg 4 limits comparable with asymmetries method 
and with the ATLAS and CMS projections. Compared 
to the method of angular asymmetries, this approach has 
an advantage of using the complete shape information of 
CP-odd observables. It is demonstrated that some of 
the observables, that do not generate significant angular 
asymmetry in presence of significant CP-mixing, can still 
provide restrictive f gi limits when their complete shape is 
analysed. Combining several CP-odd observables in the 
same fit or combining several angular asymmetries would 
likely further improve sensitivity to the CP violating cou¬ 
pling. It is noted that careful experimental investigation 
of all observables, even not the leading ones, is important, 
since they probe different terms of the HZZ vertex. 
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FIG. 7. Results of the mixing angle a scan produced with the mixing angle observable fit corresponding to the integrated 
luminosity of 300 fb -1 . 
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